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Abstract 

Context. Recent sub-millimetric observations at the Plateau de Bure interferometer evidenced a cavity at ~ 46 AU in radius into 
the proto-planetary disk around the T Tauri star LkCal5 (V1079 Tau), located in the Taurus molecular cloud. Additional Spitzer 
observations have corroborated this result possibly explained by the presence of a massive (> 5Mj up ) planetary mass, a brown dwarf 
or a low mass star companion at about 30 AU from the star. 

Aims. We used the most recent developments of high angular resolution and high contrast imaging to search directly for the existence 
of this putative companion, and to bring new constraints on its physical and orbital properties. 

Methods. The NACO adaptive optics instrument at VLT was used to observe LkCal5 using a four quadrant phase mask coronagraph 
to access small angular separations at relatively high contrast. A reference star at the same parallactic angle was carefully observed to 
optimize the quasi-static speckles subtraction (limiting our sensitivity at less than 1.0"). 

Results. Although we do not report any positive detection of a faint companion that would be responsible for the observed gap 
in LkCal5's disk (25-30 AU), our detection limits start constraining its probable mass, semi-major axis and eccentricity. Using 
evolutionary model predictions, Monte Carlo simulations exclude the presence of low eccentric companions with masses M > 6 M Jup 
and orbiting at a > 100 AU with significant level of confidence. For closer orbits, brown dwarf companions can be rejected with a 
detection probability of 90% down to 80 AU (at 80% down to 60 AU). Our detection limits do not access the star environment close 
enough to fully exclude the presence of a brown dwarf or a massive planet within the disk inner activity (i.e at less than 30 AU). Only, 
further and higher contrast observations should unveil the existence of this putative companion inside the LkCal5 disk. 

Key words. Stars: pre-main sequence - Planetary systems: proto planetary disks - Instrumentation: adaptive optics, coronagraphy 



1. introduction 

Since more than two decades, the first glimpse of planetary for- 
mation has been directly observed (Smith & Terrile 1984). Then, 
following the discovery of the first extrasolar planet around a 
solar-type star (5 1 Peg, see [Mayor & Queloz|1995[ l , more than 
400 planets have been found, mainly using the radial velocity 
technique (RV). However, the explored time span limits the ra- 
dial velocity study (as well as pulsar timing, micro-lensing, pho- 
tometric transit, astrometric techniques) to the close (< 4-5 AU) 
circumstellar environment. To understand the way exo-planetary 
systems form and evolve at wider orbits, it is then particularly in- 
teresting to use a complementary technique, such as direct imag- 
ing. 

Around young solar analogues the current deep imaging 
performances enable the detection of companions down to the 
brown dwarf (BD) and planetary mass regimes. Typical sepa- 
rations larger than 50-100 mas (i.e. > 5 - 10 AU for a star at 
100 pc) can be explored for BD companions and for planetary 
mass objects, typical separations in the range 200-400 mas (ie 
> 20 - 40 AU for a star at 100 pc). The discoveries of the plane- 



tary mass companions like 2M1207 b (Chauvin et al. 2004) and, 
more recently, around Fo malhaut (|Kalas et al.|2008"] l, a triple sys- 
tem around HR 8799 (Marois ~et al.|2009| > and a close candidate 
around /3 Pic (Lagrange et al. 2009) confirm the high contrast 
and high angular resolution imaging capabilities of current in- 
struments. To resolve ultra-cool and light companions, large tele- 
scopes equipped with adaptive optics (AO) are required, as well 
as dedicated focal devices, such as coronagraph and differential 
imagers. 

Of special interest is the detection of planets still embed- 
ded in their parental proto-planetary disk, as this may give us 
clues on the timescales for planet formation and on the disk 
dissipation mechanisms. Indirect information can be given by 
gaps and asymmetries observed in the disk structures as detected 
through scattered light images in near infrared (NIR) (e.g. for 
HD 141569, see Augere au et al.|199 9) or, through millimeter in- 
terfero metric measurements as recently done for AB Aur ( Pietu 
[etai7||2005] l and TW Hydra ( |Hughes et arT||2007> . Two contro- 
versial discoveries have been recently obtained one with RV for 
the giant planet (M ~ 10M Ju p) at 0.04 AU in the disk that sur - 
rounds TW Hya (See |Setiawan et al.|2008||Huelamo et al.|2008) , 
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the other by differential polarimetric and coronagraphic imaging 
for a point source (companion or disk over-density due to dust 
accretion on an unseen companion) in the outer annulus of dust 
of the AB Aur disk (Oppenheimer et al. 2008). More recently, 
a planetary candidate has been detected in the f3 Pic circumstel- 
lar disk (Lagrange et al. 2009). If confirmed, these companions 
will directly evidence the link between planetary formation and 
the main morphological and dynamical peculiarities of transition 
disks. 

In this paper, we present the outcomes of NACO obser- 
vations, taken using a four quadrant phase mask coronagraph, 
which aims at offering enhanced detection performances at small 
angular separations. Our goal was to detect the companion re- 
sponsible for the large gap, observed at mm wavelengths, in 
the disk that surrounds the young star LkCa 15 (see Pietu et al. 
[2006l l. 

After summarizing the properties of this star and its proto- 
planetary disk in section 2, we present the NACO observations 
and the data reduction and analysis in section 3 and 4 respec- 
tively. We finally report our detection limits in section 5 and 
discuss them in regards of the domain of mass and orbital pa- 
rameters explored for a putative companion. 



2. LkCa 15 

LkCa 15 (V1079 Tau: V = 12.09, K = 8.16, K5 and d ~ 140 pc) 
is a T Tauri star aged of ~ 3-5 Myr (see Simon et al.| 2000; 
|Bertout & Genova|2006|l, loca ted in a hole of the Taurus molec- 
ular cloud. IPietu et al.l d20061l led sub-mm observations of the 



that a ~ 5 - 10Mj up planet orbiting at 30 AU would be suffi- 
cient to evacuate the inner 50 AU of the LkCa 15 disk. Similar 



proto planetary disk of LkCa 15, using the Plateau de Bure inter- 
ferometer with an angular resolution of about 0.4". The observa- 
tions were made in "track-sharing", to observe LkCa 15 and the 
Herbig star MWC 480 with a common calibration curve. Any 
morphological difference between the two targets is thus gen- 
uine. They found a clear inner hole around LkCa 15 while MWC 
480 is centrally peaked. This is confirmed by the modeling of the 
data, with a derived inner hole of ~ 46 AU (i.e. the size of our 
own Solar System). A multi isotopes analysis of CO rotational 
lines led |Pietu et aT] ( |2007| l to derive the physical properties of 
the outer circumstellar disk surrounding the star. Contrary to the 
other system studied in a similar way, no clear vertical tempera- 
ture gradient was found in the disk structure, possibly due to the 
peculiar geometry of LkCa 15 disk. 



Pietu et al. (2006 1 discussed the possible mechanisms that 



could explain the inner hole structure (not completely empty, 
as some IR excess is present above the stellar black body. See 



Bergin et al. 2004 b. Planetary formation or the presence of a low- 



mass companion seems to be plausible since the LkCal5 disk is 
quite massive. In fact Simon et al. (2000) estimated, from the 



kinematics, a total mass for the system of Md yn = 1-00 + O.1M . 
This set the maximal mass of the putative companion to 0.2 M Q , 
since the mass of LkCa 15 could hardly be less than 0.8 M Q , due 
to its spectral type. This value is in agreeement with a recent es- 
timate of the mass of the star obt ained by|Bertout et al.|(2007| l. 
Using a new distance estimate by Bertout & Genova ( 2006), they 
confirmed LkCa 15 as a kinematic member of the Taurus asso- 
ciation seen through a hole in the molecular cloud, retrieving a 
mass value of 1.12 ± 0.08M o . 

Additional results obtained by Alexander et al. (2006) also 
suggest that other mechanisms such as photo-evaporation are un- 
likely to be effective due to the large disk density, which is sig- 
nificantly higher than the one expected if the photo-evaporation 
starts propagating beyond 20-30 AU. |Pietu et al. (2006i suggest 



conclusions were reached by Espaillat et al. (2007 1 with Spitzer 
observations and the modeling of the SED of LkCa 15. Their 
analysis suggests that a gap is present in the disk of LkCa 15, 
with an inner disk going from 0.12-0.15 AU to 4-5 AU, and an 
outer disk of inner radius 46 AU, in perfect agreement with the 
findings of Pie tu et aT] ( |2006[ ). They also concluded that plane- 
tary formation or the presence of a close-in stellar or sub-stellar 
companion are the most probable explanations for the circum- 
stellar material shape around LkCal5 



3. Observations 

3.1. Telescope and instrument 

The observations were performed on December 26th, 2007 at 
ESO/Paranal, using NACO, the AO-assisted near-IR camera 
NAOS-CONICA ( [Rousset et aLl[2003] > mounted on one of the 
Nasmyth focus of the UT4 8m-telescope. Among the numerous 
NACO observing modes ( |Lenzen et al.|2003] l, the classical and 
coronagraphic imaging modes were used. 

Coronagraphic observations were performed with the four- 
quadrant phase mask (4QPM) optimized for K s band observa- 
tions. The S13 objective (FoV of 14" X 14" and plate-scale of 
13.25 mas/pixel) was chosen for a more precise centering with 
the 4QPM and a better sampling of the PSF. The 4QPM splits 
the focal plane into four equal areas, two of which are phase- 
shifted by n. As a consequence, a destructive interference occurs 
in the relayed pupil and the on-axis starlight rejected on the edge 
of the geometric pupil is filtered with a Lyot stop, being a circu- 
lar hole 90% of the pupil size. The advantage over the classical 
Lyot mask is the possibility to access inner angular separations 
lower than 0.35 " (the smallest NACO occulting mask) at rela- 
tively large contrast (see Boccaletti et al. 2004 , 2008). 

However, a significant part of the starlight is left in the focal 
plane due to uncorrected aberrations composed of a dynamical 
halo averaging over time plus a quasi-static halo corresponding 
to optical aberrations along the optical train (from telescope to 
detector). To mitigate this problem, a coronagraphic image of 
a reference star was taken just after our science target observa- 
tions, with the same instrumental settings, to serve for speckle 
calibration in an image-subtraction process. This reference star 
(BD + 22 729, V = 11.5, K = 7.9) has similar visible and NIR 
magnitudes to ensure similar AO correction and signal-to-noise 
at the detector. Moreover, observations were performed to match 
the parallactic angle with LkCal5, i.e the instrumental pupil con- 
figuration, in order to optimize the overlap of the speckle pattern 
and diffraction spikes position in the final subtracted image of 
LkCal5. 



3.2. Observing strategy 

The coronagraphic observations were preceded by a classical 
imaging sequence that provides a photometric reference. A neu- 
tral density (1.12% transmission) and a narrow band filters at 
2.17 fim were used to serve as an image quality check and 
as photometric reference for the coronagraphic observations. 
Table [T] summarizes the observing parameters. During the coro- 
nagraphic observing sequence, the precise centering of the sci- 
ence target behind the focal plane mask was critical to max- 
imize the central star attenuation. The coronagraphic observa- 
tions were acquired at two instrument positions rotated by 33°. 



2 



M. Bonavita et al.: Searching for sub-stellar companion into the LkCal5 proto-planetary disk 



Table 1 : Observing parameters used for each source, observed 
in classical and coronagraphic imaging. 



UDject 


r vr-i 1 ^ 




DF» T~)Q 




(Classical) 


(Coronagr.) 


(Coronagr.) 


Filter 


ND + NB 2 17 


K s 


K s 


Objective 


S13 


S13 


S13 


tint (s) 


30 


24 


24 


N frame X M eX p 


2x1 


4x8 


4x8 


fi„t[sky] (s) 


30 


24 


24 


N frame X A^pfsky] 


2x 1 


4x 1 


4x 1 



Notes. The individual integration time (r int ), the number of frames 
(N frame) averaged by the detector and the number of repeated ex- 
posures (Wcxp) are reported on source and on sky. At the end, 
N f m me x ^exp gives the number of images available for the data re- 
duction and analysis. 



Sky images were immediately observed for LkCal5 and its ref- 
erence, after each coronagraphic sequence. 




4. Data reduction and analysis 

4. 1 . Image processing and selection 

The data were processed using the Eclipse [|] reduction software 
( |Devi llard 1997) for bad pixel correction, flat fielding and sky 
subtraction. Individual frames were inspected by eye to remove 
all low-quality images degraded by waffle aberrations or vari- 
able AO corrections. Finally, only LkCal5 and BD + 22 729 
images with similar parallactic angles were selected to optimize 
our PSF-subtraction. Among the initial data set of 8 images per 
source and per rotator position, only 3 were kept after the selec- 
tion process totaling an integration time of 288 seconds. 

4.2. Subtraction of the diffraction residuals 

The purpose is to properly subtract the stellar contribution from 
the LkCal5 images. For each rotator and parallactic angle po- 
sition, the reference star was shifted at a 0.1 -pixel accuracy, 
scaled and subtracted to minimize the residuals. An IDL custom- 
made tool was used to rapidly converge on an acceptable shift 
and scaling solution. Alternatively, a residual minimization us- 
ing the AMOEBA^] function was applied giving consistent re- 
sults. Fig. [T] shows the result of the single image-subtraction 
image of LkCal5 at a rotator position of 33° corresponding to 
the set that provides the best contrast. Additionally, we took ad- 
vantage of the observations taken at two rotator positions to ex- 
plore dead zones hidden by the secondary spikes and the coron- 
agraphic mask transition (see zoom-i, in Fig. 1). Derotation was 
also applied to average the speckles patterns and increase the fi- 
nal signal-to-noise. Finally, Double subtraction of images at 
and 33° were obtained to remove the non rotating aberrations 
(related to static optics in NACO) and possibly reveal a positive- 
negative signature expected for a true companion (and not for 
instrumental residuals). All sets of single, double-subtracted and 
derotated-averaged images were finally compared to derive a fi- 
nal subtracted image, optimized at different angular separations. 

1 http://www.eso.org/projects/aot/eclipse/ 



Figure 1: LkCal5 4QPM coronagraphic image at a rotator posi- 
tion of 33° after PSF subtraction and spatial filtering. The zoom- 
in non-subtracted image is reported on the upper-left corner with 
the same field of view and thresholds. 



4.3. Detection limit 

Using our final residual image, a pixel-to-pixel 2D noise map 
was estimated using a sliding box of 5 x 5 pixels over the whole 
NACO FoV (see Fig. [2] and ID-detection limits in Fig. [3]>. The 
6 cr detection limit map was obtained after renormalization by 
the LkCal5 images obtained in classical imaging and applying 
all corrections related to the use of different optical set-ups: ex- 
posure times (see Table [T}, neutral density and Lyot-stop trans- 
missions (a factor of 1.12% and 80.8% respectively, according 
to |Boccaletti et al.p 008 ) and NB2.11 to K s filter transformation. 
Finally, we also took into account that the presence of the four- 
quadrant mask causes an attenuation of the off-axis objects due 
to the 4QPM transition, decreasing our sensitivity close to the 
axis (see |Boccalett i et al.||2004) . Compared to previous deep 
imaging surveys with no a-priori information on the orbital con- 
figuration of the system, the LkCa 15 disk system properties of- 
fer the possibility to set constraints on the orbital properties of 
a possible co-planar companion. The use of the complete 2D- 
detection map is therefore mandatory to fully exploit the sky- 
projected spatial information of the LkCa 15 environment pro- 
vided by our direct imaging observations. 



5. Results 

5. 1 . Null-detection in the central hole 

At separations corresponding to the central hole (p < 0.33 ", 
46 AU) no evidence of a point-like object was found in any of our 
image-subtracted sets. The azimuthal averaged contrast places a 
limit of AK S = 7mag equivalent to about 12My for a hypothet- 
ical projected distance of 30 AU. A statistical approach of this 
detection limit is analyzed thoroughly in section 5.2 We also 



http : I /www. physics.nyu.edu/grierlab/idLhtmlJielp/A8.htmlwp992475 searched for point-like objects at separations that are not com- 
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-660. -440. -220 o. 220. 440. Figure 3: As a comparison are reported four different AK S 1D- 

RA (mas) contrast detection limits extracted at 4 position angles (45, 135, 

225 and 315 degrees). They illustrate the azimuthal variation at 
close inner angles due to the NACO PSF structure and the rele- 
Figure2: AK S contrast map giving the contrast between the vancy of using a contrast map to take that effect into account, 
faintest target detectable at a 6<x level and the primary star. 



patible with the presence of the central hole. We found a low 
significant point-source at a separation of p = (0.67 + 0.02)" at 
a position angle of PA = (340.7 + 0.3)°. Although close to the 
detection limit, this point-source lies at the boundary between 
the speckle and the background noise regimes and is then vis- 
ible in the subtracted image of Fig. [T] We coarsely estimated a 
contrast of AK S = 10.2mag in a 5 pixels aperture. Error bars are 
not estimated as the point-source is quite close to the noise level 
(less than 6<x if we refer to Fig. 3}, so the previous contrast value 
should be taken with caution. There are a number of artifacts that 
may produce such patterns like waffle and spiders. The nearest 
spider spike of which the trace is still visible in Fig.[T]is offset by 
25° while the waffle mode appears in a 45° direction at a posi- 
tion of 0.49 " therefore not compatible with the presence of this 
point-source. However, it is not visible in the subtracted image 
obtained for the 0° rotator position while at such separation the 
detection limit is almost similar for both rotator positions (see 
Fig. [2]). But the point-source would be at about 0.2 " from the 
vertical 4QPM transition and slightly attenuated. As it is difficult 
to rule out a structure in the speckle pattern, additional observa- 
tions will be mandatory to test the presence of this potentially 
low-mass object. 

5.2. Companion mass and orbital parameters 
5.2.1. Simulation description 

A first step is to convert our 6 <x detection limit map in terms 
of minimum mass map. This is usually done considering the 
star apparent magnitude in K s -bmd and its distance, finally the 
evolutionary model predictions at the age of the system to con- 
vert the absolute magnitude limits derived in the NACO filters 
in masses. We have considered here an age of 4 Myr and a dis- 
tance of 140 pc for LkCal5. We then considered two classes of 
evolutionary models based on different assumptions on the initial 
conditions: Hot Start models considering an initial spherical con- 
tracting state (Chabrier et al. 2000, Baraffe et al. 2003, Saumon 



ture of a core-accretion planet formation model ( |Marley et al. 
2007, Fortney et al.||2008] l. In the case of the core-accretion 
model predictions, our detection performances do not allow to 
access the planetary mass regime at all. Massive hot Jupiters are 
indeed predicted to be much fainter at young ages (|Marley et al. 



& Marley 2008 i and core accretion models coupling planetary 
thermal evolution to the predicted core mass and thermal struc- 



2007 ). Therefore, minimum mass maps were determined using 
Hot Stall model predictions over the planetary and brown dwarf 
mass regime. 

In a second step we used our minimum mass maps to cal- 
culate the detection probability (Pd) of companions of various 
masses and orbital parameters (semi-major axis a, eccentricities 
e, inclination i, longitude of the ascending node O, longitude 
of periastron co and time of periastron passage T p ). Considering 
the disk properties (inclination, position angle and inner radius), 
different assumptions can be made to fix partially the companion 
orbital properties. The rest of the orbital parameters can be ran- 
domly generated. Each simulated companion is then, according 
to its position on the projected orbit, placed on our 2D minimum 
mass map, and its detectability is tested, comparing its mass with 
the minimum value achievable at the same position in the FoV. 
Running the simulation 10.000 times for a given set of mass and 
orbital parameters enables to derive a detection probability. As 
already pointed out in Sec. |4.3| we took advantage of both the 
available information on the disk geometry plus the spatial in- 
formation coming from the detection limit map, in order to test 
all probable set of physical and orbital parameters for a putative 
companion. 

The smallest projected physical separation probed around 
LkCal5 is limited by the 4QPM coronagraph attenuation in- 
side 0.15" (equivalent to 3 times the angular resolution), set- 
ting the minimum semi-major axis considered in our simu- 
lations to 20 AU. Projected physical separations as large as 
1000 AU are explored, but we decided to restrain our study to 
the close circumstellar environment considering semi-major axis 
a = [20, 280] AU. Additionally, we restrain the parameters space 
explored for M - [3, 100] M Jup and eccentricity e = [0.,0.3]. 
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Figure 4: Non-detection probability map of a faint companion as a function of its mass and semi-major axis in the case of a null- 
eccentricity (circular) orbit solution. Inclination and longitude of the ascending node have been fixed using the disk properties: 
iDisk = 49° and £2 = PA Disk + 90° (= 151 ± 90). 



5.2.2. Results 

In Fig. [4] the results of our simulation are given in terms of 
non-detection probability map as a function of the companion 
mass and semi-major axis. We considered the most constrain- 
ing case of a coplanar orbit with an inclination and a longitude 
of the ascending node fixed by the disk properties: ioisk = 49° 
and Q. = PA Disk + 90° (= 151 ± 90°). Only the null-eccentricity 
map is showed, but the solutions for low-eccentricities (e < 0.3) 
are not significantly different. The range of semi-major axis es- 
timated for the predicted candidate companion (25-35 AU, see 
|Pietu et al.1|2006| l and the disk boundarie s (46 - 177 AU) are 
reported, as reported by Pietu et aL] ( |2007| l. 

Any low stellar mass companion with masses should have 
been detected with a high confidence level (non-detection prob- 
ability lower than 0.1) for semi-major axis larger than 55 AU 
(non-detection probability of 0.2 at 30 AU). Any brown dwarf 
companion with masses larger than 10 Mj up should have been 
detected for semi-major axis larger than 80 AU with a non- 
detection probability of 0.1 (60 AU for 0.2). Our non-detection 
probability map shows that we start exploring the planetary mass 
regime but not deep enough to reject the presence of any mas- 
sive giant planet at the location of the disk depletion (non detec- 
tion probability higher than 0.9 at 20-35 AU). For this physical 
separation, the presence of more massive companions, 20 Mj up 
or above, cannot be significantly excluded. At larger distances 
(a > 100 AU), the presences of companions of any mass above 
6 Mj up seem unlikely. Although we only explored low eccentric- 



ity solutions in a systematic way, the above results also provide 
lower limits on the detection probabilities for larger eccentricity 
values, as the time spent far from the star becomes progressively 
larger. 

5.2.3. Limitations 

As for all the current deep (or even wide field) imaging stud- 
ies of young, very low mass objects or companions (See e.g. 
Niels en & Close||2009t |Chauvin et aL"l|2"009"l |Lafreniere efaLl 
2007; Marois et al. 20081, the main limitations of our analysis 
come from the conversion of our detection limits in terms of 
minimum detectable masses. They are related to the uncertainty 
in the age determination of LkCal5 and the use of non-calibrated 
evolutionary models for young ages and very low masses. By 
comparison, uncertainties on the system distance and K s appar- 
ent brightness are negligible. Age and model predictions are dis- 
cussed below: 

1. LkCal5 is a confirmed member of the Taurus-Auriga asso- 
ciation, for which an age of 3-5 Myr is estimated ( |Simon| 
|et al.||2000| . To explore how the age uncertainty affects our 
results, we ran similar simulations for ages of 1, 4 and 7 Myr. 
The impact of the age uncertainty is much more important in 
the planetary mass domain than in the brown dwarf regime. 
Consequently, it does not significantly affect our conclusions 
relative to the presence of a brown dwarf or low stellar mass 
companion around LkCal5. 
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2. The applicability of evolution tracks of brown dwarfs at ages 
less than a few million years have been already cautioned by 
Baraf fe et al.| ( |2003| l. The role of the initial conditions has 
been also questioned by |Marley et al.| ( |2007| ), specifically 
for young giant planets where a connection between giant 
planet formation and evolution models seems mandatory. At 
the location of the LkCal5 inner disk cavity (between 20 and 
46 AU), our 4QPM observations are not sensitive at all to 
apparent magnitudes predicted for planetary mass compan- 
ions described by the core accretion start models. They are 
marginally sensitive for high planetary masses with fluxes 
predicted by the Hot Start models. The strongest constraints 
are actually set in the brown dwarf regime for masses larger 
than M > 20 Mj up , where the Hot Start evolutionary mod- 
els are more appropriate. Although these models still need 
to be more extensively calibrated for this range of masses 
and ages, recent discoveries of young calibrato rs fMohanty 
[eTaLl[2004l |Close et al.|[2003j |Stassun et aL|2006| l indicate 
that their predictions are relatively faithful regarding the ac- 
curacy required in our study. Our conclusions related to the 
probable presence of a brown dwarf or low mass star com- 
panion responsible for the inner hole detected in the LkCal5 
disk remain then meaningful down to 30 AU. 

6. Conclusions 

The T-Tauri star LkCa 15 was observed with VLT/NACO us- 
ing the 4QPM coronagraph, reaching a contrast lower than 9.5 
in K s band, at separations higher than 0.5 ". Our goal was the 
detection of a low-mass companion, with a mass spanning from 
0.2 M Q down to 5 Mj up , which presence has been suggested as 
an explanation for the large cavity, evidenced by sub-millimeter 
observations, in the disk surrounding the star. 

We do not report any positive detection of a close compan- 
ion to the star LkCal5. Based on our detection limits and Hot 
Start evolutionary model predictions, we ran simulations to take 
into account that the presence of a putative companion enable 
us to constrain its mass and semi-major axis given reasonable 
assumptions from the disk geometry. 

The innovative approach in our case is the combined use of 
either the available information on the disk configuration (in- 
clination, ascending node, gap prediction) and the complete 2D 
information of the coronagraphic image, to actually constrain 
the probability of existence of a stellar or sub-stellar compan- 
ion around LkCal5. As a main result of our statistical analy- 
sis we can exclude the existence of a low mass star and brown 
dwarf companion with semi-major axis larger than 55 and 80 AU 
respectively. The planetary mass regime is only partially con- 
strained and at large semi-major axis (a > 100 AU) the existence 
of massive (M > 6 M Jup ) planetary mass companions seem un- 
likely. 

The limitations of our study is also discussed in regards of 
the LkCal5 age uncertainty that does not affect much our con- 
clusions in the case of low mass star and brown dwarf com- 
panions, and of the uncertainty related to the initial conditions 
adopted for the evolutionary model, particularly critical in the 
planetary mass regime. Finally, in case of a true companion or- 
biting LkCal5 in the inner disk cavity, our observations would 
favor a planetary mass or low-mass brown dwarf companion, 
although more massive companions cannot be completely ex- 
cluded with high detection probability. 

Further deep imaging studies at the 5-50 AU scale at a new 
epoch should provide complementary information to completely 
reject the existence of close stellar or brown dwarf companion 



to LkCal5 and pursue the search for the putative companion 
that would be responsible for the disk geometry and inner cavity 
within 46 AU. 
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